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Abstract: Preliminary results of MINOS experiment indicate that #23 is not maximal. 
Global fits to world neutrino data suggest two nearly degenerate solutions for 623: one 
in the lower octant (LO: #23 < tt/4) and the other in the higher octant (HO: 623 > tt/4). 
Vn — > v e oscillations in superbeam experiments are sensitive to the octant and are capable of 
resolving this degeneracy. We study the prospects of this resolution by the current T2K and 
upcoming NOz^A experiments. Because of the hierarchy-5cp degeneracy and the octant- 
<5cp degeneracy, the impact of hierarchy on octant resolution has to be taken into account. 
As in the case of hierarchy determination, there exist favorable (unfavorable) values of 
<5cp for which octant resolution is easy (challenging). However, for octant resolution the 
unfavorable Sqp values of the neutrino data are favorable for the anti-neutrino data and 
vice-verse. This is in contrast to the case of hierarchy determination. In this paper, we 
compute the combined sensitivity of T2K and NOzaA. to resolve the octant ambiguity. If LO 
is the true octant, then NOz^A can rule out HO at 2a C.L., irrespective of the hierarchy and 
<5cp- Addition of T2K data improves the octant sensitivity. If HO and normal hierarchy 
are the true choices, then the combined data from NOz/A and T2K with its designed five 
years run in neutrino mode, is incapable of a 2a resolution of the octant for all 5qp- A 2a 
resolution of the octant, for all combinations of neutrino parameters, becomes possible if 
T2K has balanced neutrino and anti-neutrino runs of 2.5 years each. 
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1 Introduction and Motivation 

Our understanding of the smallest lepton mixing angle #13 has improved quite dramatically 
in last one year or so and finally it has been confirmed to be non-zero with unprecedented 
confidence by the reactor experiments Daya Bay [1] and RENO [2]. They have found a 
reasonably large 1-3 mixing 

sin 2 2<9i 3 | DayaB ay = 0.089 ± 0.010 (stat) ± 0.005 (syst) [3], and 
sin 2 26>i 3 | RENO = 0.113 ± 0.013 (stat) ± 0.019 (syst) [4], 

in agreement with the measurements performed earlier by T2K [5, 6], MINOS [7, 8], and 
Double Chooz [9, 10] experiments. Combined analyses of all the neutrino oscillation data 
available [11-13] imply a non-zero value of #13 at more than 10cr and predict a best-fit 
value of sin 2 #13 ~ 0.023 with a relative la precision of 10%. These recent high precision 
measurements of #13 have taken us one step further in validating the standard three flavor 
oscillation picture of neutrinos on a strong footing [14]. Also, a moderately large value of 
#13 has provided a 'golden' opportunity to directly determine the neutrino mass hierarchy 1 
(NMH) using the Earth matter effects, and to unravel the evidence of leptonic CP violation 
(LCPV) 2 in accelerator based long baseline neutrino oscillation experiments [15]. 

1 Two possibilities are there: it can be either normal (NH) if A31 = m| — ml > 0, or inverted (IH) if 
A31 < 0. 

2 If the Dirac CP phase, <5cp differs from or ty. 
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Another recent and crucial development is the indication of non-maximal #23 by the 
MINOS accelerator experiment [8]. However, the atmospheric neutrino data, dominated 
by Super-Kamiokande, still prefers the maximal value of #23 [16]. All the three global fits 
of world neutrino data [11-13] also point to the deviation from maximal mixing for #23 
i.e. |0.5 - sin 2 6*23 1 0. 

Both these new measurements, non-zero value of #13 and non-maximal #23; will provide 
crucial inputs to the theories of neutrino masses and mixings [17-20]. A number of innova- 
tive ideas, such as/iHr symmetry [21-29], A4 flavor symmetry [30-34], and quark- lepton 
complementarity [35-38] have been invoked to explain the observed pattern of one small 
and two large mixing angles in the neutrino sector. Measurements of the precise values of 
#13 and #23 will reveal the pattern of deviations from these symmetries and will lead to a 
better understanding of neutrino masses and mixings. In particular, the resolution of #23 
octant will severely constrain the patterns of symmetry breaking. 

In Ua survival probability, the dominant term is mainly sensitive to sin 2 2#23- Now, if 
sin 2 2#23 differs from 1 as indicated by the recent neutrino data, then we get two solutions 
for #23: one < vr/4, termed as lower octant (LO) and the other > 7r/4, termed as higher 
octant (HO). In other words, if the quantity (0.5 — sin 2 #23) is positive (negative) then #23 
belongs to LO (HO). This is known as the octant degeneracy of #23 [39] which is a part of 
the overall eight-fold degeneracy [40, 41], where the other two degeneracies are (#i3,<5cp) 
intrinsic degeneracy [42] and the (hierarchy, 5cp) degeneracy [43]. 

The octant ambiguity of #23 is considered to be the most difficult one to deal with 
among the eight-fold parameter degeneracies. In the past when we had only an upper 
bound on #13, a possible way of resolving this degeneracy by combining future reactor 
data with accelerator disappearance and v e appearance measurements was suggested in 
[44, 45]. Adding the information from the 'silver' channel [y e — > v T ) to the 'golden' channel 
{y e — > u^) in the proposed neutrino factory setup is demonstrated to be one of the elegant 
ways to tackle this degeneracy [46, 47]. The possibility of determining the deviation of 
#23 from maximal mixing and consequently the correct octant of #23 in very long baseline 
neutrino oscillation experiments and as well as in future atmospheric neutrino experiments 
has been discussed in [48-56]. One clear message that has been conveyed by all these 
novel works is that one can achieve a very good sensitivity to the quantity 1 0.5 — sin 2 #23) 
from the conventional beam experiments (MINOS, ICARUS and OPERA), the current 
generation superbeam experiments (presently running T2K and upcoming NOz^A), and also 
from the current (Super-Kamiokande) and future atmospheric data (India-based Neutrino 
Observatory). But, determining the sign of 1 0.5 — sin 2 #23) is deemed to be a very difficult 
job to pursue and it demands a large value of #13. 

Now, in the light of recently discovered moderately large value of #13, it would be 
quite interesting to study whether the expected appearance data from the ongoing T2K 
experiment [6, 57] in Japan and the upcoming NOi^A experiment [58, 59] in the United 
States can resolve the octant ambiguity of #23 or not? In this paper, we address this issue. 

The structure of the paper is as follows. We start in section 2 by revisiting our present 
understanding of the 2-3 mixing angle. Section 3 describes in detail the physics issues 
related to the octant of #23- We show the event rates for T2K and NOi^A in section 4. At 
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the end of this section, we also describe the simulation method followed. We present our 
results in section 5. Finally, in section 6, we summarize and draw our conclusions. Expected 
events rates in T2K and NOi^A (both for neutrino and anti-neutrino) as a function of 5cp 
can be found in Appendix A. Allowed regions in the sin 2 #23 (test) - #cp(test) plane for the 
true value of 5cp = and all the four combinations of true hierarchy and true octant are 
shown in Appendix B. 

2 Present Understanding of the 2-3 mixing angle 

Our present knowledge of 623 comes from two sources: a) atmospheric neutrinos and b) 
accelerator neutrinos. In both cases, the muon neutrino disappearance is parametrized in 
the form of two-flavor survival probability 

P w = 1 - sin 2 20 eff sin 2 ( ~ jff-J • (2-1) 

Analysis of the data gives reasonably precise values for the effective two-flavor parameters. 
Relating these to the three flavor parameters depends on the experimental set up. In the 
case of atmospheric neutrinos, the path lengths involved vary from 20 km to 13000 km and 
the energies vary from 200 MeV to a few GeV. This represents a very wide range in L/E. 
The approximations valid for some values of L/E are not valid for others. Therefore, for 
atmospheric neutrinos, it is not possible to obtain a direct relation between the effective 
two-flavor parameters and the three-flavor parameters. However, for accelerator neutrinos, 
L and E are chosen so that A.31L/E ~ tt/2. Hence, A21L/E -C 1 (A21 = — m\) and 
can be treated as a small perturbation. In this approximation, it was shown that [60-62] 

Am 2 ff = A31 — (cos 2 0i2 — cos #cp sin 0i3 sin 20i2 tan 023) A21 , (2.2) 
sin 2 20 eff = 4 cos 2 6 13 sin 2 23 (l - cos 2 6 13 sin 2 23 ) . (2.3) 

The mixing angles and phases are defined according to the Particle Data Group convention 
[63]. The atmospheric neutrino data, dominated by Super-Kamiokande, still prefers the 
maximal value of sin 2 20 e g- = 1 (> 0.94 (90% C.L.)) [16]. But the preliminary results from 
the MINOS accelerator experiment favor a non-maximal value of sin 2 20 c g = 0.94jlQQ 5 [8]. 



Reference 


Forero et.al. [11] 


Fogli et.al. [12] 


Gonzalez-Garcia et.al. [13] 


sin 2 23 (NH) 
3cr range 


n4 r, 7 +o.034 ffinfil 0+0.022 

U - 4/ ' -0.027 a? U - Did -0.040 

0.36 -> 0.68 


0.386«;° 2 t 
0.331 -> 0.637 


n 41 +0.037 ffi n trq+0.021 
U - 4i -0.025 U - t,y -0.022 

0.34 -> 0.67 


sin 2 6» 23 (IH) 
3cr range 


600 +U U26 
u.ouu_ .o 31 

0.37^ 0.67 


n 009+0.039 
u.o»z_ 022 

0.335 -> 0.663 



Table 1: la bounds on sin 2 #23 from the global fits performed in References [11], [12], and [13]. NH 
and IH stand for normal and inverted hierarchies. The numbers cited from Ref. [13] are those obtained by 
keeping the reactor fluxes free in the fit and also including the short baseline reactor data with L < 100 m, 
with the mass hierarchy marginalized. 
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Global fits, using three-flavor oscillations, give information directly on 623 rather than 
# c fj. The best-fit values and ranges of #23 from the three recent global fits [11], [12], and 
[13] are listed in Table 1. A common feature that has emerged from all the three global fits 
of the world neutrino data is that we now have indication for non-maximal 623- Thus, we 
have the two degenerate solutions: either 623 belongs to the LO (sin 2 623 ~ 0.4) or it lies 
in the HO (sin 2 623 ~ 0.6). This degeneracy, in principle, can be broken with the help of 
V/m ^ v e oscillation data. The preferred value would depend on the choice of the neutrino 
mass hierarchy. However, as can be seen from Table 1, the fits of reference [11] do not 
agree on which value should be preferred, even when the mass hierarchy is fixed to be NH. 
In [12], LO is preferred over HO for both NH and IH. Reference [13] marginalizes over the 
mass hierarchy, so the degeneracy remains. The global best-fits in references [11, 13] do 
not observe any sensitivity to the octant of Q23 unless they add the atmospheric neutrino 
data. But, in [12], they do find a preference for LO even without adding the atmospheric 
data. In this paper, we take the best-fit value of sin 2 #23 in the lower octant (LO) to be 
0.41 while that in the higher octant (HO) to be 0.59 [13]. 

3 Physics of the octant of Q 2 3 

In a long baseline experiment, za, charged current (CC) events are the most copious. These 
experiments can measure j/„ — > fn survival probability, as a function of energy. The 
reconstruction of the minimum of leads to precise values for |Am 2 ff | and sin 2 26 c g 
[8, 60-62]. Therefore, we get two degenerate best fit values for sin 2 823, one in LO and the 
other in HO, with small allowed regions around them. The P w expression has subleading 
terms which are sensitive to octant [64]. But, these are suppressed by the small parameter 
a = A21/A31 and their biggest impact occurs for energies where P w is very small. So, the 
overall octant sensitivity of P w is negligible. However, P w gives a precise measurement 
of sin 2 29 e g which in turn gives precise allowed regions for sin 2 #23 m the two octants. 

In the presence of matter, the —> v e oscillation probability, expanded perturbatively 
in a, can be written as [64-66] 

P^e = ft sin 2 023 + ft cos(A + 5 CP ) + ft cos 2 623- (3.1) 

In the above equation, P^ e is written in a way to highlight the octant and <5cp dependent 
terms. Here 



I- A) 2 

sinAi sinA(l - A) 



ft = a cos 6<i3 sin 2^12 sin 2^i3 sin 2^23 



A I - A 



/33 = a 2 sin 2 2^ 1 2Cos 2 1 3 Sin ^ 4 , (3-2) 

A 2 

with A = A31L/4S, A = A/A31. A is the Wolfenstein matter term [67] and is given 
by ^4(eV 2 ) = 0.76 x 10~ 4 /? (g/cc)£?(GeV). p is the density of matter in the Earth. For 
NOi^A and T2K, this is set equal to the density in the crust of 2.8 g/cc. 
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For normal hierarchy (NH), A31 is positive and for inverted hierarchy (IH), it is nega- 
tive. The matter term A is positive for neutrinos and is negative for anti-neutrinos. Hence, 
in neutrino oscillation probability, A is positive for NH and is negative for IH; vice- verse for 
anti- neutrinos. Moreover, sign of 5qp is reversed for anti-neutrinos. The left(right) panel 
of figure 1 shows .vs. Fi u (Pp S .vs. Ep) for all possible values of 5qp and for the two 
values of sin 2 #23 , assuming NH to be the true hierarchy. These plots are for the experiment 
NOz/A . 




E [GeV] E [GeV] 

Figure 1: P Me as a function of neutrino energy. The left panel (right panel) is for v (y). Here, the bands 
correspond to different values of Sep from — n to tt. These plots are for NOi/A (L=810 km), sin 2 20 13 = 0.089 
and NH. For LO(HO), sin 2 6 2 z = 0.41(0.59). 

As can be seen from these figures, there is an overlap between the two octant bands 
for some values of 5qp- This overlap spoils the octant determination sensitivity of the 
experiment. Those combinations of octant and <5cp which lie farthest from the overlap 
will be the favorable combinations while the overlapping combinations will be the unfa- 
vorable ones. For favorable combinations, octant degeneracy can easily be resolved while 
it is difficult for unfavorable ones. A comparison of left and right panels leads to a very 
interesting observation. For the v case, the favorable combinations are (HO, 5qp ~ —90°) 
and (LO, 5qp ~ 90°). However, for the v case, the the favorable <5cp values for HO are 
~ 90° while those for LO are ~ —90°. Similar arguments hold for IH also. This is in 
contrast to the hierarchy-<5cp degeneracy, where for a given hierarchy, the favorable <5cp 
region is the same for both v and v. Thus, we are led to the conclusion that a balanced 
neutrino and anti- neutrino data is imperative for resolving the octant for all values of <5cp- 

Let us do a small, quantitative analysis of the octant-<5cp degeneracy of P^ e . For 
simplicity, here we keep the hierarchy fixed. Increase in #23 increases P^ e . While a change 
in <5cp can increase or decrease P^. For different <5^p and ^p, P Ate (LO, d}jp) may be very 
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close to P (Je (HO, Sqp ). For the degenerate case, P IMe (LO,6}yp) = P Me (HO, Sq®), leading to 
cos(A + <j£g) - cos(A + 5%$) = ^ ~ - (sin 2 - sin 2 0$). (3.3) 

P2 

For the NOi^A baseline L = 810 km and energy of peak flux E = 2 GeV, we get for NH 
and v 

cos(A + 6%$) - cos(A + = 1.7. (3.4) 
The above equation will have solutions only if 

0.7 < cos(A + <f£g) < 1-0, 
-1.0 < cos(A + 5^§) < -0.7. (3.5) 

From this, we get their ranges to be: 

- 116° < 6$ < -26°, 

64° < 5^ < 161°. (3.6) 

Thus, we find that for NH and v of energy 2 GeV, P Me (LO,-116° < 5 C p < -26°) is 
degenerate with P^HO, 64° < 5c~p < 161°). A similar analysis can be done for T2K 
baseline of L = 295 km and energy of peak flux E = 0.8 GeV. The overlap regions in 
the ^cp range are essentially the same as above because the values of A and those of 
{fix — fiz)/ $2 are nearly the same for the two experiments. This is to be contrasted with 
the hierarchy discrimination, where the overlap range for NOz^A is very different from that 
of T2K because of widely different matter effects [68, 69]. 

In figure 2, we show the values of <5cp> which lie in the overlapping region, as a function 
of energy, for the combinations LO-NH and HO-NH. For NOi/A experiment, most of the 
signal events come from the range 1.5 - 2.5 GeV. Hence, we consider only this energy 
range in the figure. The blue-dotted (red-crossed) shaded region shows those values of 5qp, 
for which there is an overlap between Pj^P and Pj^P, for v (p) data. These plots show 
the degenerate octant-^cp values in the relevant energy range and help explain the $cp 
dependence seen in octant sensitivity. 

The values of <5cp which lie farthest from any of the overlap regions will be the most 
favored ones. Those <5cp values which lie in one of the overlap region (i.e. either in v or 
i>) but are far from the other overlap region will also be favored (though to a lesser extent 
than the previous values) if balanced v and v runs are taken. The regions of <5cp which 
are common or very close to both overlap regions will be the most unfavored ones. We 
find that 5cp ~ is the most unfavorable for LO-NH and HO-IH whereas <5cp ~ 180° is 
the most unfavorable for HO-NH and LO-IH. This pattern is observed in the results in 
subsection 5.1. 
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Figure 2: Plots showing octant-<5cp degeneracy in P Me as a function of neutrino energy. The left (right) 
panel is for LO (HO). The blue-dotted (red-crossed) regions are for v (v). For a given E^E^), P Me [ LO, 
vertical blue-dotted (red-crossed) <5cp region in the left panel ] values are degenerate with P^e\ HO, vertical 
blue-dotted (red-crossed) b~cp region in the right panel] values. The exact degenerate octant-<5cp values can 
be found out using equation 3.4. As an example, for E„ of 2 GeV, P Me (LO, <5cp = —90°) is degenerate 
with P Me (HO, S C p « 66°). These plots are for NO^A (L=810 km), sin 2 28 13 = 0.089 and NH. For LO(HO), 
sin 2 023 = 0.41(0.59). 



4 Event Rates for T2K and NOz/A 



In this paper, we simulate the data for the two off-axis superbeam experiments: T2K 
and NO^A using GLoBES [70, 71]. In the T2K experiment, a v a beam from J-PARC is 
directed towards Super-Kamiokande detector, 295 km away. The flux peaks sharply at 0.8 
GeV, close to the first oscillation maximum in P{y a — > u e ). The experiment is scheduled 
to run for 5 years in the neutrino mode only. The details of T2K experiment are given in 
[57]. The information regarding signal efficiencies and backgrounds are taken from [72, 73]. 
NOz^A is a 14 kT totally active scintillator detector located at Ash River; a distance of 810 
km from Fermilab. The flux peaks at 2 GeV, again close to the first oscillation maximum 
in P{v a v e)- This experiment is scheduled to have three years run in neutrino mode first 
and then later, three years run in anti-neutrino mode as well. The details of the experiment 
are given in [74]. In light of the recent measurement of large #13, NOi/A has reoptimized 
their signal and background acceptances. In our calculations, we use these reoptimized 
values, the details of which are given in [59, 75]. 

In our simulations, we have used the following input values for neutrino oscillation 
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parameters [8, 13]. 



|Am 2 ff | = 2.4 x l(T 3 eV 2 

A 21 = 7.5 x 10" 5 eV 2 , sin 2 12 = 0.3 

sin 2 26»i3 = 0.089 

The value of A31 is calculated separately for NH and for IH using equation 2.2 where 
Am 2 ff is taken to be +ve for NH and -ve for IH. The uncertainties in the above parameters 
are taken to be o"(|Am 2 ff |) = 4% [8] and cr(sin 2 26*13) = 5% [76]. The solar parameters (A21 
and sin 2 O12) and the Earth matter density in the calculation of matter term are held fixed 
throughout the calculation. 




Figure 3: v and v appearance events for all possible combinations of hierarchy, octant and Sep- The left 
(right) panel is for T2K (NOz/A). Here sin 2 29 w = 0.089. For LO(HO), sin 2 6 23 = 0.41(0.59). Note that for 
T2K, equal v and v runs of 2.5 years each has been assumed. 



In figure 3 , we plot v events .vs. v events for various octant-hierarchy combinations. 
In each case, with varying values of <5cp> the plot becomes an ellipse. The left panel shows 
these ellipses for T2K whereas the right panel shows the same for NO/vA . Here, we assumed 
that T2K will have equal v and v runs of 2.5 years each. In the right panel, we see that 
the ellipses for the two hierarchies overlap whereas the ellipses of LO are well separated 
from those of HO. Hence, we can expect that NOz^A will have better octant discrimination 
capability than hierarchy discrimination. This situation is even more dramatic in the left 
panel where there is large overlap between the two hierarchies but clear separation between 
the octants. Thus, it is very likely that v data from T2K may play an important role in 
the determination of octant. 

3 The variation of v and v appearance events with <5cp is shown in appendix A for both T2K and NOz^A . 
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We illustrate the octant determination capability of NOzM. by considering some special 
points in the right panel. The features that we emphasize here can also be discerned from 
figure 2. 

• LO-NH, <5cp = 0°: The coordinates of this point in the v-v event plane are (70,23). 
No point on the HO-NH ellipse has these coordinates but there exist a number of 
points with coordinates close to these. Therefore, the wrong octant is difficult to rule 
out for this point. 

• LO-NH, #cp = 90°: The coordinates of this point in the v-v event plane are (50,26). 
The v events are much lower than those of all points on the HO-NH ellipse but there 
is degeneracy in the v events. Hence, ruling out the wrong octant should be possible 
for this point. 

• LO-NH, 5cp = 180°: The coordinates of this point in the v-v event plane are (62,18). 
Both v and v events are much lower than the corresponding events of any point on 
the HO-NH ellipse. Thus, ruling out the wrong octant will be the easiest for this 
point. 

• LO-NH, <5cp = —90°: An argument similar to the case 5cp = 90° can be made, except 
that the v events have degeneracy between octants, but the v events for this point 
are much below those of any point on the HO-NH ellipse. 

Similar arguments can be made for other octant-hierarchy combinations with the ex- 
ception that the most favorable and most unfavorable 5qp values will differ. 

We see that for the given choice of best-fit values of sin 2 623 , NOz^A has very good 
sensitivity to octant resolution due to balanced v and v runs. Because the favorable and 
unfavorable values of Sqp (pertaining to octant resolution) are different for v and v, no 
5cp value is absolutely unfavorable if balanced v and v runs are taken. Therefore, a 
single experiment on its own, can have very good sensitivity as illustrated in figure 3. 
This is in stark contrast to the case of hierarchy, where we saw that having data from 
two experiments with widely different baseline is a necessity because the favorable and 
unfavorable combinations turned out to be the same for v and v. 

Before discussing our results, we briefly describe the numerical procedure adopted. We 
calculate Ax 2 using the default definition in GLoBES which is Poissonian. We minimize 
this Ax 2 to compute the octant resolution capability. For a given octant-hierarchy com- 
bination and a true value of 5cp, we compute the events spectra and label it data. Then 
we compute the theoretical events spectra where the octant is chosen to be the wrong 
one and the neutrino parameters are randomly chosen within their allowed 3a ranges. We 
then calculate the Ax 2 between the data and the theoretical spectra. This calculation uses 
the Ax 2 defined in GLoBES, which is valid for a Poissonian distribution. We add a Ax 2 
coming from gaussian priors on | Am 2 ff | and on sin 2 2^13. The systematic uncertainties are 
included using the method of pulls. If the Ax 2 ^ > 4, then we can say that the wrong 
octant is ruled out at 2a for the given combination and the given true value of 5cp- This 
calculation is repeated for all true values of 5qp and for all combinations. If the Ax 2 ^ > 4 
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in each case then the wrong octant can be ruled out independently of true Sqp, hierarchy 
and octant. 



5 Results 



In this section, through various plots, we show the sensitivity of T2K and NOi^A to the 
octant of #23- 

5.1 Allowed regions in test sin 2 823 - true Sqp plane 

In figures 4-7, we have plotted the values of test sin 2 623 allowed by T2K and NOzM. data 
as a function of true <5cp for each of the four combinations of octant and hierarchy. In our 
calculations, we do not assume a prior knowledge of hierarchy and hence consider both the 
possibilities for test hierarchy. Octant can be determined only if the wrong octant values 
are ruled out for both possibilities. 

In this analysis, we have included the data from both the disappearance channel P w 
and the appearance channel P^e. We varied sin 2 623 within its 3a range: [0.34, 0.67] allowed 
by the current global fits. Here, we are constraining sin 2 623 only and are marginalizing 
over 5qp and hence take the 2a limit relevant for 1 d.o.f. The contours in these figures are 
defined by Ax 2 < 4. The disappearance channel gives a precise measurement of sin 2 29 e g 
which leads to two narrow, <5cp -independent, allowed sin 2 623 bands, one in each octant. 
The appearance channel, because of its large octant sensitivity, discriminates against the 
wrong octant band. If the statistics are large enough, the wrong octant can be ruled out. 



LO-NH true, NH test, 2a 



LO-NH true, IH test, 2a 



NOvA 
NOvA+T2K[5+0] 




0.32 




0.32 




5 CP (true) 



8 rp (true) 



Figure 4: Allowed values of test sin 2 #23 at 2a (1 d.o.f.) C.L. as a function of true Sep. LO-NH is 
assumed to be the true octant-hierarchy combination. The left (right) panel corresponds to NH (IH) being 
the test hierarchy. 
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LO-IH true, NH test, 2a 



LO-IH true, IH test, 2a 




5 rp (true) 



8 rp (true) 



Figure 5: Allowed values of test sin 2 #23 at 2a (1 d.o.f.) C.L. as a function of true Sep- LO-IH is assumed 
to be the true octant-hierarchy combination. The left (right) panel corresponds to NH (IH) being the test 
hierarchy. 



HO-NH true, NH test, 2a 



HO-NH true, IH test, 2a 




5 CP (true) 



8 rp (true) 



Figure 6: Allowed values of test sin 2 823 at 2a (1 d.o.f. ) C.L. as a function of true Sep. HO-NH is 
assumed to be the true octant-hierarchy combination. The left (right) panel corresponds to NH (IH) being 
the test hierarchy. 
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Figure 7: Allowed values of test sin 2 623 at 2a (1 d.o.f. ) C.L. as a function of true S C p. HO-IH is 
assumed to be the true octant-hierarchy combination. The left (right) panel corresponds to NH (IH) being 
the test hierarchy. 
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Figure 8: Allowed values of test sin 2 823 at 2a (1 d.o.f. ) C.L. as a function of true Sep- HO-NH is 
assumed to be the true octant-hierarchy combination. The left (right) panel corresponds to NH (IH) being 
the test hierarchy. Note that for T2K, equal v and v runs of 2.5 years each has been assumed. 



From figures 4 and 5, we note that the combined data of T2K and NOi^A can rule out 
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HO if LO is the true octant. This holds for both NH and IH. However, if HO is the true 
octant, then LO can not be ruled out for a reasonably large fraction of true 5qp- 

As argued in section 3, the favorable and unfavorable regions are very different for 
neutrinos and anti-neutrinos. Hence, we explored if an improvement in the octant deter- 
mination can be achieved if T2K has equal neutrino and anti- neutrino runs of 2.5 years 
each 1 . As expected, in the two cases of LO-NH and LO-IH, a (2.5^ + 2.5z7) runs in T2K 
is also effective in ruling out the wrong octant at 2a. Hence, we have not displayed the 
corresponding figures. But, for the two cases HO-NH and HO-IH, the balanced v and v 
runs are more effective than the (5f) run. This is shown in figures 8 and 9. In the left 
panel of figure 8, there is a very small sliver of allowed sin 2 #23 in the wrong octant for true 
<5cp in the range [130°, 160°]. But, the A X 2 for this sliver is very close to 4. Hence, this 
wrong octant region can be effectively discriminated against, leading to a determination of 
the octant of 623 for any true value of Sep ■ 

The above set of figures also give the precision on sin 2 623 that one can obtain. This 
precision is the result of the precise measurement of sin 2 2#23 from the disappearance chan- 
nel. From figures 4 and 5, we get <5(sin 2 #23) = 0.015 if LO is the true octant. From figures 
6 and 7, we see that a 5 year v run of T2K combined with NOz^A gives <5(sin 2 <9 2 3) = 0.02 
if HO is the true octant. When the T2K run is changed to (2.5z/ + 2.5P), the precision in 
sin 2 #23 becomes slightly worse because of the loss of statistics caused by v run as seen in 
figures 8 and 9. 



HO-IH true, NH test, 2a 



HO-IH true, IH test, 2a 




5 rp (true) 



8 CP (true) 



Figure 9: Allowed values of test sin 2 823 at 2<r (1 d.o.f. ) C.L. as a function of true Sep- HO-IH is 
assumed to be the true octant-hierarchy combination. The left (right) panel corresponds to NH (IH) being 
the test hierarchy. Note that for T2K, equal v and v runs of 2.5 years each has been assumed. 

4 In GLoBES, the background events for T2K are given for a 5 year v run taken from [73]. We have 
taken care to do appropriate scaling of these events in computing for (2.hv + 2.5z>) runs. 
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5.2 Ax 2 .vs. true 5qp plots 



In this section, we study the behavior of A% 2 between the true and the wrong octants as a 
function of true <5cp- Here, the A% 2 is computed in the following way. First, we fix the true 
value of <5cp- We take sin 2 #23 to be its best-fit value in the true octant: 0.41 for LO and 
0.59 for HO. If the LO (HO) is the true octant, the test values of sin 2 9 23 in the HO (LO) 
are varied within the range [0.5, 0.63] ([0.36, 0.5]), where 0.63 (0.36) is the 2a upper (lower) 
limit of the allowed range of sin 2 823- The A% 2 is computed between the spectra with the 
best-fit sin 2 823 of the true octant and that with various test values in the wrong octant 
and is marginalized over other neutrino parameters, especially the hierarchy, sin 2 2#i3 and 
5qp- Figures 10 and 11 show the minimum of this A% 2 .vs. the true value of <5cp. 

From figure 10, we see that the NOi^A data by itself can almost rule out the wrong 
octant at 2a, if LO is the true octant. We see that the A% 2 dips just below 4 for true 
Sqp ~ 0(180°) if the true hierarchy is NH (IH). But, as argued earlier, this small allowed 
region can be effectively discriminated because of the relatively large Ax 2 . If HO is the 
true octant, then NOi^A data is not sufficient to rule out the wrong octant as seen in figure 
11. In fact, the wrong octant can be ruled out only for about half of the true 5cp values. As 
illustrated in figures 10 and 11, addition of T2K data improves the octant determination 
ability significantly. 




Figure 10: Octant resolving capability as a function of true <5cp for various set-ups . In these plots, LO 
is assumed to be the true octant. The left (right) panel corresponds to NH (IH) being the true hierarchy. 
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Figure 11: Octant resolving capability as a function of true <5cp for various set-ups . In these plots, HO 
is assumed to be the true octant. The left (right) panel corresponds to NH (IH) being the true hierarchy. 



From figure 10, we see that the combined data from NOi^A and T2K (5v) give a 2<t 
octant discrimination for all values of true 5cp if LO is the true octant. From figure 11, 
we see that this combined data can rule out the wrong octant at 2a for HO-IH, but not 
for HO-NH. The problem of HO-NH can be solved if the T2K has equal v and v runs of 
2.5 years each. This change improves the octant determination for the unfavorable values 
of true <5cp (where A% 2 is minimum) for all four combinations of hierarchy and octant. In 
particular, for the case of HO-NH, it leads to a complete ruling out of the wrong octant at 
2o" for all values of true 5cp. Thus, balanced runs of T2K in v — v mode is preferred over 
a pure v run because of better octant determination capability. 

Figures 10 and 11 show that the combined data from NOi^A and T2K has a better 
overall octant resolving capability if LO is the true octant. We found out that this feature 
of LO being more favorable compared to HO is a consequence of marginalization over the 
oscillation parameters (mainly 5cp) an d the systematic uncertainties. We checked that in 
the absence of any kind of marginalization A%^ Q is consistently larger than A^lq • 

Before we summarize, we would like to address a non-trivial issue. As mentioned 
earlier, MINOS favors a non-maximal value of sin 2 2t? e ff = 0.94. #23 values calculated from 
this sin 2 26 c g value will not be complimentary and both LO and HO values will be closer 
to 45°, compared to what we have in [13]. Octant resolving capability with these inputs 
for #23 will be more challenging. This issue is under study [77]. 
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6 Summary and Conclusions 



Recently, the preliminary results from MINOS experiment have indicated that #23 is not 
maximal. This raises the question of the true octant of 623 i.e. whether 623 < vr/4 (LO) 
or 623 > vr/4 (HO). v e appearance searches at the presently running T2K and upcoming 
NOz^A experiments are sensitive to the octant, especially in light of moderately large #13. 
The main difficulty in #23-octant resolution stems from octant-5cp degeneracy. We explored 
this degeneracy in detail and found that balanced neutrino and anti-neutrino exposures 
are mandatory to overcome this problem. This is because the octant-(5cp combinations, 
which have degenerate probabilities in the neutrino data are free of this degeneracy in the 
anti-neutrino data and vice-verse. We also studied the prospects of T2K and NOzaA. to 
improve the precision in 6*23. The disappearance data leads to two degenerate solutions 
of sin 2 #23, one in each octant. If LO (HO) is the true octant, the expected precision will 
be <5(sin 2 #23) = 0.015 (0.02). The v e appearance data can resolve the octant degeneracy. 
NOz^A alone can rule out the wrong octant at 2a if LO is true, independently of hierarchy 
and 5cp- Combined data from T2K and NOz/A can establish the correct octant at 2a for 
any combination of hierarchy, octant and 5qp provided T2K has equal 2.5 years of exposure 
in both neutrino and anti-neutrino modes rather than a full 5 years neutrino run. 
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A Events .vs. 5 C p 



In this appendix, we consider the variation of appearance events as a function of 5qp in 
both v and v modes. 
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Figure 12: Total appearance events rates for all possible combinations of octant and hierarchy as a 
function of the Sep- The left (right) panel is for v (p) running. These plots are for NOz^A (L=810 km), 
sin 2 2(9i3 = 0.089. For LO(HO), sin 2 6> 23 = 0.41(0.59). 




Figure 13: Total appearance events rates for all possible combinations of octant and hierarchy as a 
function of the <5cp. The left (right) panel is for v (u) running. These plots are for T2K (L=295 km), 
sin 2 26>i3 = 0.089. For LO(HO), sin 2 6> 23 = 0.41(0.59). 

In figure 12, we show the variation of Un — > v e appearance events .vs. <5cp for NOzaA. in 
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both v and v modes. We see from the left panel {y events) that the combinations HO-NH 
and LO-IH are well separated but the other two combinations- HO-IH and LO-NH have 
essentially the same event numbers. But in the right panel (p events), HO-IH and LO-NH 
are well separated and the other two combinations are nearly degenerate. Thus, we see 
that the unfavorable combinations in v mode are favorable in v mode and vice-verse. This 
feature is seen in the event rates for T2K also as shown in figure 13 where we have assumed 
T2K to have 2.5 years of each neutrino and anti- neutrino run . 



B Allowed regions in test <5cp - test sin 2 6*23 plane 



In section 3, we saw that the CP conserving values of <5cp pose the biggest challenge in 
octant determination. Therefore, in this appendix, we simulate T2K and NO^A data 
for true 5cp = and analyze it. Figures 14-17 are drawn for each of the four possible 
combinations of true hierarchy and true octant. These figures show the regions allowed by 
the data in the test <5cp - test sin 2 823 plane at 2a (Ax 2 < 6.18 for 2 d.o.f.). 
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Figure 14: Allowed regions in test <5cp - test sin 2 623 plane at 2a (2 d.o.f.) C.L. for true Sep = 0. 
LO-NH is assumed to be the true combination. The left (right) panel corresponds to NH (IH) being the 
test hierarchy. 
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Figure 15: Allowed regions in test Sep - test sin 2 823 plane at 2a (2 d.o.f.) C.L. for true Sep = 0. 
LO-IH is assumed to be the true combination. The left (right) panel corresponds to NH (IH) being the test 
hierarchy. 
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Figure 16: Allowed regions in test Sep - test sin 2 823 plane at 2a (2 d.o.f.) C.L. for true Sep = 0. 
HO-NH is assumed to be the true combination. The left (right) panel corresponds to NH (IH) being the 
test hierarchy. 
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Figure 17: Allowed regions in test Sep - test sin 2 823 plane at 2a (2 d.o.f.) C.L. for true Sep = 0. 
HO-IH is assumed to be the true combination. The left (right) panel corresponds to NH (IH) being the 
test hierarchy. 

These figures illustrate that, in general, the combined data of T2K and NOz^A discrim- 
inate against the wrong octant of #23- Also, note that, these figures give us the allowed 
ranges in both sin 2 623 and 5qp (2 d.o.f.) around the right and wrong octant for the single 
value of true 5qp = 0. 
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